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A Method and Apparatus Based on Combination of Physiological Parameters 
for Assessment of Analgesia during Anesthesia or Sedation 

5 FIELD OF THE INVENTION 

The present invention relates to a method and an apparatus for monitoring a condi- 
tion of a patient under anesthesia or sedation, while, in particular, information ac- 
quired from different sources in the patient is utilized. 

10 

BACKGROUND OF THE INVENTION 

Concept of the depth of anesthesia has been of interest for recent decades, and sev- 
eral measures have been proposed to assess the depth of anesthesia. Recently, how- 

1 5 ever, this unitary anesthesia theory of the existence of one-dimensional concept 
called "depth of anesthesia" has been strongly criticized as oversimplified. Instead it 
has been suggested that the anesthesia has not one but three main components: hyp- 
nosis, analgesia and muscle relaxation. Different anesthetic regimens have different 
effect on these three components. Furthermore, they have effects on both cortical 

20 and sub-cortical levels. An adequate anesthesia means unresponsiveness to both 
noxious and non-noxious stimuli. The former may be defined by means of hemody- 
namic, motor and endocrine stability, while the latter is related to the loss of con- 
sciousness and recall and amnesia. In practice the adequate anesthesia is adminis- 
tered by using a combination of drugs with different effects on brain, spinal cord, 

25 autonomic nervous system and neuromuscular junction. The combination of these 
effects hence creates the hypnotic, analgesic and muscle relaxing effects. 

In general anesthesia the patient is conducted through the phases of anesthesia from 
the induction to the varying lengths of maintenance period and to the final emer- 

30 gence out from anesthesia. Though the patient does not usually recall any surgical 
events or perceive surgical pain, the recovery and the post-operative comfort of the 
patient very often depend on the quality of the anesthesia during the operation itself. 
Adequate administration of analgesic drugs - meaning that over-doses and under- 
doses can be avoided during anesthesia - is believed to advance the recovery of the 

35 patient. It has been suggested that this is due to two main reasons. Surgical pain 
may sensitize the pain pathways during surgery and thus lower the pain threshold in 
such a way that even rather intense pain management in the post-operative period is 
ineffective. It is said that the best way to avoid post-operative pain is a good and 
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adequate administration of analgesics during operation. The other mechanism is 
probably related to the secretions of stress hormones during surgery. These hor- 
mones may have their effects long after surgery and can slow down the physical and 
psycho-physiological heeling of the patient. Adequate pain medications can sup- 
5 press the autonomic nervous system and prevent excess secretion of these stress 
hormones. In this context the term "nociception" is commonly used to refer to the 
perception of pain. The receptors involved in pain detection are aptly enough re- 
ferred to as nociceptors. Nociceptive input is conveyed from the peripheral end or- 
gans to the central nervous system. Projection neurons in the spinal dorsal horn pro- 

10 ject to cell nuclei in supraspinal areas such as the thalamus, the brainstem, the mid- 
brain etc. Of these, the synaptic junctions in the thalamus play a very important role 
in the integration and modulation of spinal nociceptive and non-nociceptive inputs. 
Nociceptive inputs are finally conducted to the cortex, where the sensation of pain 
is perceived. Stimulation of these central nervous system regions either electrically 

15 or chemically, e.g. by morphine and other opiates, produces analgesia in humans. 

Currently the anesthesia practices rely on rather subjective assessments of the ade- 
quacy of the drug treatment during anesthesia. Anesthesiologists observe the patient 
and decide for the proper drugs they give to the patient. Though this often is enough 

20 to avoid adverse events such as arousal or muscle movements during surgery, which 
in fact very seldom occur in normal anesthesia nowadays, more objective measures 
for the anesthesia are needed. Recently the progress in the biopotential signal analy- 
sis has lead in reasonable quantitative estimation of the hypnotic level of the patient, 
and thereby the titration of the anesthetic agents can be guided by these new meas- 

25 urements. 

For instance, the neurological activity of the brain is reflected in biopotentials avail- 
able on the surface of the brain and on the scalp. Thus, efforts to quantify the extent 
of anesthesia-induced hypnosis have turned to a study of these biopotentials. The 

30 biopotential electrical signals are usually obtained by a pair, or plurality of pairs, of 
electrodes placed on the patients scalp at locations designated by a recognized pro- 
tocol and a set, or a plurality of sets or channels, of electrical signals are obtained 
from the electrodes. These signals are amplified and filtered. The recorded signals 
comprise an electroencephalogram or EEG, which normally has no obvious repeti- 

35 tive patterns, contrary to other biopotential signals, like electrocardiogram (ECG). 
Among the purposes of filtering is to remove electromyographic (EMG) signals 
from the EEG signal. EMG signals result from muscle activity of the patient and 
will appear in electroencephalographic electrodes applied to the forehead or scalp of 
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the patient. They are usually considered artifacts with respect to the EEG signals. 
Since EMG signals characteristically have most of their energy in a frequency range 
from 40 Hz to 300 Hz, which is different than that of the EEG, major portions of the 
EMG signals can be separated from the contaminated EEG signal. 

5 

A macro characteristic of EEG signal patterns is the existence of broadly defined 
low frequency rhythms or waves occurring in certain frequency bands. Four such 
bands are recognized: Alpha waves are found during periods of wakefulness and 
may disappear entirely during sleep. The higher frequency Beta waves are recorded 
10 during periods of intense activation of the central nervous system. The lower fre- 
quency Theta and Delta waves reflect drowsiness and periods of deep sleep. 

For clinical use, it is desirable to simplify the results of EEG signal analysis of the 
foregoing, and other types, into a workable parameter that can be used by an anes- 

1 5 thesiologist in a clinical setting when attending the patient. Various such parameters 
for relating EEG signal data to the hypnotic state of the patient are discussed in the 
literature. Several use frequency domain power spectral analysis. These parameters 
include peak power frequency (PPF), median power frequency (MPF), and spectral 
edge frequency (SEF). A peak power frequency (PPF) parameter uses the frequency 

20 in a spectrum at which occurs the highest power in the sampled data as an indication 
of the depth of anesthesia. The median power frequency (MPF) parameter, as its 
name implies, uses the frequency that bisects the spectrum. In the same fashion, the 
spectral edge frequency uses the highest frequency in the EEG signal. To improve 
the consistency of an indicator of the hypnotic state or depth of anesthesia, several 

25 parameters are often employed in combination. For example, the spectral edge fre- 
quency (SEF) parameter may be combined with the time-domain burst suppression 
ratio (BSR) parameter to improve the consistency and accuracy with which the 
depth of anesthesia can be indicated. Also more complex combinations of parame- 
ters, like bispectral index (BIS) have been described. 

30 

There are a number of concepts and analytical techniques directed to the complex 
nature of random and unpredictable signals, like EEG. One such concept is entropy. 
Entropy, as a physical concept, describes the state of disorder of a physical system. 
Applying the concept of entropy to the brain, the premise is that when a person is 
35 awake, the mind is full of activity and hence the state of the brain is more nonlinear, 
complex, and noise like. Since EEG signals reflect the underlying state of brain ac- 
tivity, this is reflected in relatively more "randomness" or "complexity" in the EEG 
signal data, or, conversely, in a low level of "order." As a person falls asleep or is 
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anesthetized, the brain function begins to lessen and becomes more orderly and 
regular. As the activity state of the brain changes, this is reflected in the EEG sig- 
nals by a relative lowering of the "randomness" or "complexity" of the EEG signal 
data or conversely, increasing "order" in the signal data. When a person is awake, 
5 the EEG data signals will have higher entropy and when the person is asleep the 
EEG signal data will have a lower entropy. 

According to the International Publication WO-02/32305 both the EEG and EMG 
signal data are typically obtained from the same set of electrodes applied, for exam- 

10 pie, to the forehead of the patient. The EEG signal component dominates the lower 
frequencies (up to about 30 Hz) contained in the biopotentials existing in the elec- 
trodes and EMG signal component dominates the higher frequencies (about 50 Hz 
and above). The presence of EMG signal can provide a rapid indication of the con- 
scious-unconscious state of the patient. Importantly, because of the higher fre- 

1 5 quency of the EMG data signal, the sampling time can be significantly shorter than 
that required for the lower frequency EEG signal data. This allows the EMG data to 
be computed more frequently so that the overall diagnostic indicator can quickly in- 
dicate changes in the state of the patient. In one embodiment of the publication, the 
EEG signal data and the EMG signal data are separately analyzed and thereafter 

20 combined into a diagnostic index or indicator. As noted above, because of the celer- 
ity with which changes in the anesthetic state of the patient can be determined from 
the EMG data, the overall index can quickly inform the anesthesiologist of changes 
in the state of the patient. 

25 Another way is to observe a photoplethysmographic (= PPG) signal, which is ob- 
tained by measuring the intensity of light transmitted through or reflected by the tis- 
sue. The dynamic part of the signal is caused by variations in blood volume and per- 
fusion of the tissue, affecting scattering and absorption of the incident light. The 
most usual application of the signal is the measurement of the oxygen saturation of 

30 blood. The pulse waveform of the PPG signal is closely similar to that of the intra- 
arterial blood pressure. The waveform is reflecting the interaction between left ven- 
tricular output, i.e. cardiac output or stroke volume, and the capacitance of the vas- 
cular tree, also called vascular resistance. Blood pressure is determined by the car- 
diac output, which is stroke volume multiplied by heart rate, and vascular resis- 

35 tance. However, in addition to these global circulatory parameters, the dynamic ca- 
pacitance of the vasculature affects also the nonlinear relationship of PPG signal 
and circulatory parameters. Especially complex is the relationship between the PPG 
waveform shape within one pulse and the integrated pulse-to-pulse variables. The 
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PPG signal is related to the changes in peripheral tissue blood volume and blood ab- 
sorptivity. As it is the blood flow, that causes the blood volume changes, the PPG 
signal is hence indirectly related to local blood flow. The flow, in turn, depends on 
the pressure gradient and local vascular dynamic resistance and capacitance. 

5 

The PPG measuring as such has been utilized for a long time. For instance United 
States Patent 6,117,075 discloses a method and device for monitoring the depth of 
anesthesia (= DOA) during surgery by analyzing patterns and characteristics of os- 
cillatory phenomena in measured pulse pressure and skin temperature signals. The 

10 method utilizes pulse pressure and skin temperature oscillatory patterns describe the 
nature of sympathetic vasomotor tone. The method monitors DOA in two ways. 
Spectral characteristics of skin temperature or pulse pressure oscillatory phenomena 
are used to describe the depth of anesthesia, and the concordance between oscilla- 
tory patterns of two physiological signals, which have been recorded simultaneously 

15 but at different locations, are used to describe the depth of anesthesia. According to 
the publication a PPG signal of an anesthetized patient is continuously monitored, 
and the recorded raw PPG signal is then processed so as to generate a signal depict- 
ing the beat-to-beat pulse pressure amplitude. Then the signal is derived by detect- 
ing peaks, and calculating the difference between each positive-negative peak pair, 

20 after which a further signal is processed in a manner so as to derive a data set de- \ 
scribing very low frequency variations in pulse pressure over time in the 0.01-0.04 
Hz range, that is, the PPG signal amplitude variability. Power spectrum analysis is 
finally performed on said further signal, and the received frequency power spectrum 
characteristics are used to describe the DOA, such that a progressively narrower 

25 bandwidth describe a progressively deeper level of anesthesia. 

However, it appears that the position of the dicrotic notch as well as the PPG ampli- 
tude are dependent on various other sources than the status of vasoconstriction or 
vasodilatation, including fluid balance, temperature of the site of PPG measurement, 
30 heart rate, etc. Hence, these parameters may be interpreted with caution. Further- 
more, they refer to the usage of the PPG information as a measure of the depth of 
anesthesia, which is an oversimplified one-dimensional assumption as described 
above. 

35 The publication E. Seitsonen, M. van Gils, I. Korhonen, K, Korttila, A. Yli-Hankala: 
"EEG, Heart Rate, Pulse Plethysmography and Movement Responses to Skin Inci- 
sion" - A-582, 2002 ASA Meeting Abstracts, October 16, 2002, discloses studies 
concerning evaluation of analgesia and nociception during general anesthesia. Raw 
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EEG, the bispectral index, electrocardiography (ECG) and PPG data were collected 
and analyzed offline. RR-interval (RRI) tachogram was derived from ECG and 
frontal electromyography (fEMG) from EEG, and several beat-to-beat morphology 
parameters were derived from PPG signal, as well as various time and frequency 
5 domain parameters were computed from EEG, RRI and PPG data. When derived 
variables calculated as ratios or differences between post-incision and pre-incision 
values were compared, RRI, amplitude of the dicrotic notch in PPG, EEG spectral 
entropy and fEMG power appeared as primary variables in the optimal linear dis- 
criminant function between movers and non-movers. Finally it was concluded that 
10 combination of these parameters may be useful in assessing the level of analgesia 
and nociception during anesthesia. However, the publication does not provide any 
practical procedure for monitoring a patient. 

United States Patent 6,338,713, discloses a system and method for providing infor- 
mation to the user of a medical monitoring or diagnostic device to aid in the clinical 
decision making process. The preferred embodiment uses two estimators or predic- 
tors of the same physiological quantity, with each of the estimators being designed 
to detect specific states or artifacts in the estimated parameter and thus operating at 
a different point on its respective ROC curve; one chosen to provide high sensitiv- 
ity, the other chosen to provide high specificity. The divergence between the esti- \ 
mators is indicated by the use of a shaded region between their respective time 
trends. The use of two estimators of the same parameters with different performance 
characteristics allows the system and method of the present invention to derive ad- 
ditional information about the underlying physiologic process over and above that 
which would be available from a single estimator. The system and method of the 
present invention can derive information from not only the instantaneous values of 
the estimators and the difference between them, but also from the time trend of the 
difference. Accordingly, this publication is primarily directed to the accuracy of a 
measurement, but does not discuss the problem how the level of analgesia could be 
reliably monitored while the patient is under anesthesia or sedation. 

Anyway, the administration of analgetics is still largely based on the visual observa- 
tions of the vital signs and the hemodynamic responses of the patient to surgical 
stimulation. Analgesic drugs are usually given, when the heart rate or blood pres- 
35 sure show fast increases or are in long term at the high end of the normal ranges. 
Different motoric responses, sweating and lacrimation of the patient can be ob- 
served as well. A further problem not considered in the publications is the possible 
suppressive effect of relaxants on at least some signals acquired from patient, which 
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can have adverse effect on the reliability of the results received. The concept of an- 
algesia is very complex and due to inter- and intraindividual variability the com- 
bined specificity and sensitivity for the single-parameter based methods is not very 
good. 

5 

Accordingly the main object of the invention is to achieve a method and apparatus 
for monitoring the anesthesia or sedation of a patient so that a reliable data about 
level or depth of analgesia would be available to an anesthetist or to other purposes. 

10 The second object of the invention is to achieve a method and apparatus for moni- 
toring the anesthesia or sedation capable of using measured signals derived from 
various sources of the patient, which means that the method should not be depend- 
ent on any single type of detector. 

15 The third object of the invention is to achieve a method and apparatus for monitor- 
ing the anesthesia or sedation capable to deliver such results as an output, with the 
basis of which the adequacy of analgesia could be reliably enough assessed by in- 
experienced anesthetists or other operators, too. 

20 The fourth object of the invention is to achieve a method and apparatus for monitor- 
ing the anesthesia or sedation functioning with an acceptable speed so that a change 
in analgesia to a hazardous direction is detected and reported early enough to allow 
timely corrective actions. 

25 SUMMARY OF THE INVENTION 

According to the first aspect of the invention the method for monitoring a condition 
of a patient under anesthesia or sedation, comprises the steps of: acquiring in real- 
time at least a first signal representing a cardiovascular activity of the patient; deriv- 

30 ing, from said first signal, continuously at least a first and a second instantaneous 
parameter value related to a quantity selected from a group of quantities including 
waveform amplitudes, waveform periodicity, waveform morphology, and waveform 
variability; applying a predetermined mathematical index for probability of patient 
comfort, in which function said at least first and second parameters are variables; 

35 calculating successively changing values of said mathematical index; and indicating 
said successive index values. 
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In this case the first signal concerning cardiovascular activity is measured non- 
invasively using preferably photoplethysmography, though a pressure metering can 
be also used, whereupon the quantity for the first parameter value is a pulse wave 
amplitude, or a dicrotic notch height in the pulse wave, and the quantity for the sec- 
5 ond parameter value is a pulse rate, or a heart beat interval, or a temporal position of 
the dicrotic notch. 

As an alternative, the method according to the invention comprises the steps of: ac- 
quiring in real-time at least a first signal and a second signal representing a cardio- 

10 vascular activity of the patient; deriving, from said first signal and second signal, 
continuously at least a first and a third instantaneous parameter value related to a 
quantity selected from a group of quantities including waveform amplitudes, wave- 
form periodicity, waveform morphology, and waveform variability; applying a pre- 
determined mathematical index for probability of patient comfort, in which function 

15 said at least first and third parameters are variables; calculating successively chang- 
ing values of said mathematical index; and indicating said successive index values. 

In this case said first signal concerning cardiovascular activity is also preferably 
measured non-invasively using photoplethysmography, or measured using a pres- 
20 sure metering, whereupon the quantity for said first parameter value is a pulse wave 
amplitude, or a dicrotic notch height in the pulse wave. The second signal concern- 
ing cardiovascular activity is a cardiac excitation measured non-invasively using 
electrocardiogram, whereupon the quantity for said third parameter value is a heart 
rate of the electrical excitation. 

25 

As a further alternative, the method according to the invention comprises the steps 
of: acquiring in real-time at least a first signal and a third signal representing a car- 
diovascular and respectively a combined electrical biopotential on skull activity of 
the patient; deriving, from said first signal, continuously at least a first instantane- 

30 ous parameter value related to a quantity selected from a group of quantities includ- 
ing waveform amplitude, waveform periodicity, waveform morphology, and wave- 
form variability; calculating, from said third signal waveforms, continuously at least 
a fourth instantaneous parameter value related to a quantity selected from a group of 
quantities including energy, power, signal complexity and frequency content, each 

35 over a predetermined time period; applying a predetermined mathematical function 
for probability index of patient comfort, in which function said at least first and 
fourth parameters are variables; calculating successively changing values of said 
mathematical index; and indicating said successive index values. 
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In this case the first signal concerning cardiovascular activity is measured non- 
invasively using photoplethysmography, just as above, and the quantity for said first 
parameter value is a pulse wave amplitude, or a dicrotic notch height in the pulse 
5 wave. The third signal is a non-invasive measurement concerning neuromuscular 
and brain activity comprising electromyography and electroencephalogram. The 
electric myographic component is then extracted from said third signal, which 
originally is a combination of EMG and EEG, as an EMG partial signal and electro- 
encephalography component as an EEG partial signal. Then the quantity for said 
10 fourth parameter value is a spectral power calculated from said EMG partial signal, 
and the quantity for said fourth parameter value is a subtraction of a responsie en- 
tropy, calculated from said third signal as a whole, and a state entropy, calculated 
from said EEG partial signal. 

15 The first signal concerning cardiovascular activity can also be a cardiac excitation 
measured non-invasively using electrocardiogram, whereupon the quantity for the 
first parameter value is a heart rate or inter beat interval of the electrical excitation. 

In the context of this alternative it is also possible, and prefarable as is now be- 
20 lieved, to acquiring online a second signal representing a cardiovascular activity of 
the patient; deriving, from said and second signal, a third instantaneous parameter 
value related to a quantity selected from a group of quantities including waveform 
amplitudes, waveform periodicity, waveform morphology, and waveform variabil- 
ity; and introducing said at least third parameter in said predetermined mathematical 
25 function as an additional variable. Here the quantity for said first parameter value is 
a heart rate or inter beat interval of the electrical excitation. 

In all cases the first and the second parameter values, or the first and the third pa- 
rameter values, or the first and the fourth parameter values respectively are normal- 

30 ized on the basis their respective parameter values acquired continuously over a 
predetermined fixed time window including or excluding the latest real-time pa- 
rameter value. The above mentioned parameter values for normalizing are acquired 
either from said patient prior to incision, or prior to intubation, or prior to starting 
anesthesia or sedation, or from a group of patients prior to or during incision and/or 

35 intubation and/or anesthesia or sedation. The mathematical index for probability is a 
nonlinear equation, or a neural network algorithm, or a defined or fuzzy rule-based 
reasoning procedure. 
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The successively calculated value or values is/are indicated or informed or made 
available to the operating person or persons for further actions, and/or possibly in- 
dicated or informed or made available to an additional apparatus for further process- 
ing. It shall be noted that the absolute value of the probability index according to the 
5 invention is the most valid data made available to the personnel, though a change of 
the index may increase informativeness, too. 

According to the second aspect of the invention the apparatus for monitoring a con- 
dition of a patient under anesthesia or sedation, comprises at least first sensor means 

10 for online receiving substantially continuous electrical signal representing a cardio- 
vascular activity of the patient; first time measuring means and a first voltage/cur- 
rent dependent circuit connected with said sensor means; first memory means stor- 
ing criteria of predetermined pulse wave parameters to be extracted from said sig- 
nal; first deriving means connected to said first memory means, said first time 

15 measuring means and said first voltage/current dependent circuit for continuously 
extracting first and second values of said predefined pulse wave parameters; second 
calculation means for successively performing a predetermined mathematical pro- 
gram having temporally variable probability index values, based on said first and 
second pulse wave parameter values, as an output; and a display and/or connections 

20 into further devices. 

As an alternative, the apparatus according to the invention comprises at least first 
and second sensor means for online receiving substantially continuous electrical 
signals representing a cardiovascular activity of the patient; first and second time 

25 measuring means and a first and second voltage/current dependent circuit connected 
with said sensor means; first memory means storing criteria of predetermined pulse 
wave parameters to be extracted from said signals; first deriving means connected 
to said first memory means, said first time measuring means and said first volt- 
age/current dependent circuit for continuously extracting first values of said prede- 

30 fined pulse wave parameters; second deriving means connected to said first memory 
means, said second time measuring means and said second voltage/current depend- 
ent circuit for continuously extracting third values of said predefined pulse wave pa- 
rameters; second calculation means for successively performing a predetermined 
mathematical program having temporally variable probability index values, based 

35 on said first and third pulse wave parameter values, as an output; and a display 
and/or connections into further devices. 
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As a further alternative, the method according to the invention comprises at least 
first and third sensor means for online receiving substantially continuous electrical 
signals representing a cardiovascular activity and respectively an electrical biopo- 
tential activity of the patient; first and third time measuring means and a first and 
5 third voltage/current dependent circuit connected with said sensor means; first 
memory means storing criteria of predetermined parameters to be extracted from 
said signals; first deriving means connected to said first memory means, said first 
time measuring means and said first voltage/current dependent circuit for continu- 
ously extracting first values of predetermined pulse wave parameters; first calcula- 

10 tion means connected to said first memory means, said third time measuring means 
and said third voltage/current dependent circuit for continuously extracting fourth 
values of predetermined biopotential parameters; second calculation means for suc- 
cessively performing a predetermined mathematical program having temporally 
variable probability index values, based on said first and fourth parameter values, as 

15 an output; and a display and/or connections into further devices. 

Generally speaking, at least one height or amplitude parameter, i.e. parameter in y- 
direction in an orthogonal coordination, is derived or calculated from the signal or 
signals for usage in the function of probability according to the invention, and at 
20 least one time or temporal parameter, i.e. parameter in x-direction in an orthogonal 
coordination, is also derived or calculated from the signal or signals for usage in the 
function of probability according to the invention. In this context the entropy pa- 
rameter is considered to a height or amplitude parameter. 

25 The main advantage of the method as compared to the state-of-art is that by the 
method it is possible to significantly increase the specificity of analgesia monitoring 
as compared to a usage of a single variable as an input. The method is also practical 
in the sense that it is based on the new usage of the signals which are either rou- 
tinely monitored during anesthesia, or which may be easily monitored with the ex- 

30 isting patient monitors, e.g PPG, EEG etc., and hence the method does not require 
expensive new sensors or monitoring means to be attached to the patient. 

By combining information from various sources each related to level of analgesia 
but alone providing insufficient sensitivity and/or specificity for analgesia assess- 
35 ment in practice, it is possible to improve the quality of analgesia level monitoring. 
The most potential input signals comprise heart rate (as measured from ECG, pulse 
plethysmography, blood pressure or some other signal related to functioning of 
heart and providing beat-to-beat rhythm of the heart), plethysmography pulse pa- 
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rameters (for example plethysmography pulse notch position information) and fron- 
tal electromyography signal (fEMG). 

BRIEF DESCRIPTION OF THE DRAWINGS 

5 

FIG. 1 shows an example of a typical non-invasive PPG signal, or non-invasive or 
invasive blood pressure metering signal, representing cardiovascular pulse wave or 
blood pressure of a patient, which signal can be used as a starting raw signal in the 
invention. 

10 

FIG. 2 shows an example of a typical non-invasive ECG signal representing cardiac 
excitation of a patient, which signal can be used as a starting raw signal in the in- 
vention. 

15 FIG. 3 shows an example of a typical non-invasive EEG+EMG signal representing 
a combination of brain activity and neuromuscular activity of a patient, i.e. biopo- 
tentials, which signal can be used as a starting raw signal in the invention. 

FIG. 4 shows a typical power spectrum of a biopotential signal EEG+EMG. 

20 

FIGS. 5A - 5D are an example of real data recorded from a patient representing the 
intermediate results and final probability index calculated according to the inven- 
tion. Fig. 5D provides annotations related to surgical procedures. 

25 FIGS. 6-8 show the main steps of the preferred first, second and third embodiment, 
respectively, of the invention for processing the signal acquired from a patient. 

FIGS. 9-11 show the main components of the preferred first, second and third em- 
bodiment, respectively, of the invention as flow chart. The components shown can 
30 be separate electronic units, but quite as well portions of a computer program or 
programs for microprocessors. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

35 In the first embodiment of the invention the combination of the parameter values 
used in the probability function P for the patient comfort in anesthesia are obtained 
only from one cardiovascular signal SI. Though the outmost simplicity of this 
method in practice, multiple of quantities can be derived in the way that they still 



13 



carry different type of information from the function of the autonomic nervous sys- 
tem (= ANS). The heart rate and its variability, which mainly reflects the parasym- 
pathetic activity, can be extracted from the periodicity of the measured pulse wave, 
i.e. from the beat-to-beat pulse rate. On the other hand the amplitude of the pulse 
5 wave, especially when measured in peripheral tissue, carries information from the 
skin circulation, which is mainly controlled by the sympathetically branch of ANS. 
In case when the parasympathetical activity of the patient is stronger than the sym- 
pathetical activity, the heart rate is low, blood pressure is normal and peripheral cir- 
culation works normally as the blood vessels are dilated. If the patient has pain or is 

10 uncomfortable, the heart rate is higher and blood vessels are more constricted, 
which results in smaller peripheral pulse wave amplitudes, especially in the plethys- 
mography wave. These signs of excess sympathetical activity result also in elevated 
blood pressures easily measured by continuous blood pressure meters. In summary, 
the sympathetical activation, which increases heart rate, induces vasoconstriction in 

15 the peripheral circulation and increases blood pressure, is an indication of patient 
discomfort. Usually this appears as a sudden change in the plethysmography 
waveform parameters and blood pressure. On the other hand the low activity of the 
sympathetical and parasympathetical branch of ANS and the balance between these 
regulation mechanisms is an indicator of adequate patient comfort. Accordingly, it 

20 can be understood that even one signal can carry quite a detailed information from 
the state of the autonomic nervous system and from the patient comfort. 

In the simplest embodiment only oxygen saturation measurement with the pulse 
wave signal (PPG) is enough to calculate the probability index of patient comfort. 
25 As saturation measurement is in practice compulsory in anesthesia this advanta- 
geous embodiment of the invention is almost always available to all patients. 

In the second embodiment of the invention the combination of the parameter values 
used in the probability function P for the patient comfort in anesthesia are obtained 

30 only from two cardiovascular signals SI and S2. The heart rate is extracted from the 
ECG waveform instead of the pulse wave of the photoplethysmographic or blood 
pressure measurement. This is advantageous, because the heart beat interval can be 
determined more accurately, and more detailed information can be extracted espe- 
cially from the parasympathetical activity and from the parasympathetical — sympa- 

35 thetical balance of ANS. The better precision of the beat-to-beat interval is due to 
typically higher sampling frequencies of the ECG measurement and more impor- 
tantly due to the fact that the R-peak in the ECG waveform is very much sharper 
than the pulse wave peak. When the signal quality in PPG, due to poor blood circu- 
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lation, is bad or when the PPG or BP waves are extremely round, the exact beat in- 
terval can still be determined with precision from the ECG waveform. As the preci- 
sion of the determination of the R-R interval from ECG — even without poor patient 
condition - is always very good, more parameters, e.g. heart rate variability (HRV), 
5 can be extracted from the temporal cardiac function. These improve the perform- 
ance of the probability index for the patient comfort and leads into better sensitivity 
and specificity of the index. 

In the third embodiment of the invention the combination of the parameter values 

10 used in the probability function P for the patient comfort in anesthesia are obtained 
only from a cardiovascular signal SI and a biopontential signal S3 on skull or scalp. 
This is one of the preferred embodiments, together with the embodiment utilizing 
three signals SI, S2 and S3, regarding the specificity and sensitivity of combining 
parameters into one index of patient comfort, is achieved with the biopotential 

15 measurement on patient skull. The cardiovascular measurement are good indicators, 
but sometimes too unspecific to patient comfort as they are influenced by other psy- 
chological emotional triggers or by physiologic factors such as the intra-vascular 
volume status of the patient, by sepsis, arrhythmias, respiratory disturbance or hy- 
pertension. Therefore, in addition to the cardiovascular quantities, information is 

20 needed from the brain activity, its psychophysiological aspects and from the overall \ 
state of the central nervous system. By biopotential measurement on skull, including 
electroencephalography, electromyography and skin conductivity measurement, one 
can derive information, which is specific to the cortical brain activity and hypnosis, 
to the facial muscle nerve activity and to the sudomotor (sweating) activity, respec- 

25 tively. Hypnosis describes the state of human cortex, especially various states of 
consciousness and unconsciousness, i.e. ability of cortical processing of the sensory 
nerve information. Facial muscle nerve activity describes the motoric control of 
these muscles (from cortex), but also unintentional facial expressions, which origi- 
nate from the brain stem level and reflect also the state of the autonomic nervous 

30 system or the state of the central nervous system at the subcortical level. Finally the 
sudomotor activity is emotionally triggered from the cortex, but as the facial 
'mimic' muscles, its activation is often a result of painful or discomfortable stimula- 
tion without real perception of pain at the cortical level. It is well known that dis- 
comfort increases sweating - during anesthesia it is a common practice to sense the 

35 patient forehead for sweating, which among anesthesiologist is often interpreted as 
a sign of discomfort even in deep hypnosis. It is believed that both the sudomotor 
and electromyographic activity are specific indicators of the sympathetical activa- 
tion and patient discomfort and that the cortical activity is an additional parameter 
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for the overall suppression in the central nervous system. All this information is 
complementary within itself, but in combination with the cardiovascular signals, it 
forms an entity, by which the patient discomfort and comfort can be estimated even 
in deep hypnosis. Even more than that, the index of hypnosis and the index of anal- 
5 gesia can be made rather orthogonal to each other. This means that the indices are 
independent and that the anesthesiologist can make specific decision about whether 
patient needs more analgesic or hypnotic drugs. 

If the biopotential measurement is arranged on the forehead of the patient, both the 
10 facial muscle (fEMG) and brain (EEG) activity components are present in the raw 
measured signal S3. In principle, the same arrangement is suitable to measure the 
skin contact of the electrodes with the scalp, in which the sweating of the forehead 
is as well seen. We, however, omit this contribution and concentrate on the signal 
with EEG and fEMG, only. We also limit our discussion to one of the advantageous 
15 embodiment of the invention, namely to a technique of spectral signal entropy, the 
concept of which can be applied to both EEG and fEMG signal components. 

Entropy, when considered as a physical concept, is proportional to the logarithm of 
the number of microstates available to a thermodynamic system, and is thus related 

20 to the amount of "disorder" in the system. In this context, entropy describes the ir- 
regularity, complexity, or unpredictability characteristics of a signal. In a simple ex- 
ample, a signal in which sequential values are alternately of one fixed magnitude 
and then of another fixed magnitude has an entropy value of zero, i.e. the signal is 
completely regular and totally predictable. A signal in which sequential values are 

25 generated by a random number generator has greater complexity and higher en- 
tropy. Entropy is an intuitive parameter in the sense that one can visually distin- 
guish a regular signal from an irregular one. Entropy also has the property that it is 
independent of absolute scales such as the amplitude or the frequency of the signal: 
a simple sine wave is perfectly regular whether it is fast or slow. In an EEG applica- 

30 tion, this is a significant property, as it is well known that there are interindividual 
variations in the absolute frequencies of the EEG rhythms. 

There are various ways to compute the entropy of a signal. In frequency domain, 
spectral entropy may be computed in ways generally known. The starting point of 
35 the computations is the spectrum of the signal. Such an algorithm is implemented in 
the Datex-Ohmeda Entropy® Module, in which the discrete Fourier transformation 
of the signal is calculated. A flat spectrum indicates a high and peaky spectrum a 
low entropy value. Without going into more details, we next describe how the 
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physical concept of spectral entropy is employed in clinical decision making about 
the patient comfort. 

Referring to fig. 4, it is informative to consider two entropy indicators, one over the 
5 EEG dominant frequency range alone and another over the complete range of fre- 
quencies, including both EEG and EMG components. The State Entropy (SE) is 
computed over the frequency range from 0.8 Hz to 32 Hz. It includes the EEG- 
dominant part of the spectrum, and therefore primarily reflects the cortical state of 
the patient. The Response Entropy (RE) is computed over a frequency range from 

10 0.8 Hz to 47 Hz. It includes both the EEG-dominant and EMG-dominant part of the 
spectrum. It is advantageous to normalize these two entropy parameters in such a 
way that RE becomes equal to SE when the EMG power, i.e. sum of spectral power 
between 32 Hz and 47 Hz, is equal to zero, as the RE-SE-difference then serves as 
an indicator for EMG activation. Consequently, RE varies from 0 to 1, whereas SE 

15 is always smaller than 1 . When there is EMG activity, spectral components within 
the range 32-47 Hz differ significantly from zero and RE is larger than SE. With 
these definitions, SE and RE both serve their own informative purposes for the an- 
esthesiologists. The State Entropy is a quantity which is stable enough to provide 
the anesthesiologist at one glance of the number an idea of the current cortical state 

20 the patient is in. The Response Entropy, on the other hand, reacts fast to changes as 
it can be calculated using shorter time windows than SE. 

There are several parameters that could be classified as describing waveform mor- 
phology, including at least: the area under the pulse, FWHM (= full width at half 

25 maximum value of the pulse), amplitude/period , dicrotic notch height in absolute 
units, dicrotic notch height in relative units (= height/amplitude), skewness of the 
pulse, any characteristics after derivation or multiple of derivations of the signal, 
any parameter after FFT analysis (= first harmonic peak/base frequency peak, etc), 
any signal processing or modelling parameter (= e.g. AutoRegressive-model of the 

30 signal). In this context the parameters that could be classified as describing wave- 
form amplitude include: the relative and absolute difference between the peak value 
and the minimum value of the pulse wave. It is especially pointed out that here the 
"wave" includes a bias, generally called as a dc-component of the signal, which 
variates considerably slower than the pulse, i.e. the ac-component of the signal. Ac- 

35 cordingly, there is pulsation of the signal and a level of the signal, around which 
level said pulsation occurs. There are some parameters that could be classified as 
describing waveform periodicity, including at least: pulse rate, heart rate, heart beat 
interval, temporal position of the dicrotic notch in absolute units (after the peak- 
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position or after any fixed position in one single pulse) or in relative units (= tempo- 
ral position/period), P-R interval, Q-T interval, QRS duration. There are a plurality 
of parameters that could be classified as describing waveform variability (= changes 
of any waveform characteristics with time), including at least: SD (== standard de- 
5 viation over a fixed time window), SSD (= SD of the successive signal points= SD 
after derivation of the signal), RMS-value (as above), LF (== low frequency variabil- 
ity of any waveform characteristics: is related with the sympathetical activity of 
ANS), HF(= high frequency variability of any waveform characteristics: is related 
with the parasympathetical activity of ANS), Entropy of any characteristics (meas- 

10 ure of disorder, see Spectral entropy of EEG). It is also noted here that non-station- 
arity of the signal is one cause for waveform variability and the indices and meth- 
ods, by which the non-stationarities are detected, produce variables, which shall be 
interpreted as waveform variability. Non-stationarity here is any statistically signifi- 
cant temporal deviation of statistical measures of the signal from the average values. 

15 As to the signal complexity, there are at least the following parameters: Spectral en- 
tropy, approximate entropy, Lempel-Ziv Complexity, Golmogorov- Sinai entropy, 
Fractal exponents, bispectral indices, correlation lengths (in time). At least the fol- 
lowing parameters can be classified in frequency content: FFT spectrum derived pa- 
rameters like spectral entropy, power at different brain rhythm bands, median fre- 

20 quency, spectral edge frequency, 95% power edge frequency, mean frequency, peak 
amplitude frequency (= highest power frequency). The power and energy are self- 
explanatory parameters, and does not need further explanation. 

In the figures 1 to 3 the following parameters, which are the most frequently used 
25 and preferred parameters, are visualized: pulse wave amplitude y A ; dicrotic notch 
height yN in the pulse wave; heart beat interval r A , temporal position of the dicrotic 
notch r N ; systolic blood pressure F s ; and diastolic blood pressure 7 D . The entropies, 
which are frequently used and preferred parameters, too, cannot be visualized. The 
parameters listed above are all familiar to a person skilled in the art so that it is not 
30 necessary to mark them in the figures. 

The normalizing discussed above means that the prevailing or present or latest pa- 
rameter value is subtracted from the calculated average value of the historical pa- 
rameter values, or vice versa, or the prevailing or present or latest parameter value 
35 is divided by the calculated average value of the historical parameter values, or vice 
versa. One of these, two of these, or three of these, or all of these can be used simul- 
taneously in the predetermined mathematical function according to the invention. 
The average value can be a constant average value over time period, or over a num- 
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ber of heart beat pulses, or a moving weighted or non- weighted average value over 
a time period or over a number of heart beat pulses during operation of the patient, 
or a constant group average over certain patient type(s). 

5 The mathematical index for probability can be generally formed, according to the 
invention, several ways. Also the history data, i.e. the parameter values acquired 
continuously over a predetermined fixed time window including or excluding the 
latest real-time or prevailing parameter value, can be processed several ways. The 
probability index of the invention combines in real-time the historical parameter 

10 values with the prevailing parameter data from at least two sources. The index P or 
P(t) for probability of patient comfort is a function, which shall be understood to be 
either an equation, or an algorithm, or a reasoning or logical procedure, in which at 
least two parameters are variables This mathematical combination may be done by 
using some nonlinear equation, logical rules and/or operators, artificial neural net- 

15 works, fuzzy logic, etc. 

The index of patient comfort may be calculated by the following nonlinear equa- 
tion: 

20 P(t) = - - ^ (aX(t)+bnt)+cZ{n+d) 
where 

P is the index for probability of patient comfort, 

a, b, c, d are parameters selected to optimise the performance of the index, 

25 X(t), Y(t) t Z(t) are monitored current values of the normalized parameter values, 
t is the sample time. 

It should be noted that the coefficients of the mathematical equation must be se- 
lected according to the quantities and parameters used and according to the avail- 
30 able signal, and will depend on the exact choice and implementation of the moni- 
tored ensamble of parameter history values and their normalization means. 

The index of patient comfort may be also determined by using rule-based reasoning. 
These rules are by translating clinical knowledge acquired in our experimental stud- 
35 ies on nociception into formal "if... then... else" statements that can be directly im- 
plemented into a computer programming language. The probability index P is then 
calculated as a sequence of reasoning statements. It is possible to use exactly de- 
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fined rules as well as fuzzy rules in the reasoning. The value ranges of the parame- 
ters as well as the index for probability are in this case discredited into certain inter- 
vals to allow for symbolic reasoning. For example the index P can be subdivided 
into three or more intervals; if 0 < P < 0.3 -> P is labeled as having "low" value, if 
5 0.3 < P < 0.7 -» P is labeled as having "medium" value, and if 0.7 < P < 1 .0 -> P is 
labeled as having "high" value. Similar divisions could be done for the monitored 
values X(t), Y(t), Z(t). For example: 

if (<X(0 - "high") or (X(t) = "medium")) then 
P(t) := "low" 

else 

if ((7(0 = "low") and (Z(0 - "high")) then 
P(t) := "high" 

else 

P(t) := "medium" 

where 

= comparison operator (test of equality), 

:= assignment operator (the value of the operand on the right hand side is as- 
20 signed to the operand on the left hand side). 

This example indicates that when X(i) has a medium or high value then the prob- 
ability of nociception P is low. lfX(t) has a low value and 7(0 is low and Z(t) has a 
high value then P is high. In other cases P is medium. 

Further, the index of patient comfort may be also determined by using previously 
recorded data that has been labeled at each time instance with an indication of 
whether there was nociception or not, e.g., annotated by an anesthetist. This kind of 
procedure can be called neural network. Typical examples cases of X(t), Y(t), and 
Z(t) as recorded during instances of nociception and non-nociception can then be 
constructed. This can be done by grouping vectors v(0 = ( X(t), Y(t) t Z(t) ) that have 
been measured during nociception periods together into a set of i "typical" nocicep- 
tion examples, {vz} , and grouping vectors of data that have been measured during 
non-nociception periods into a set of j typical non-nociception examples, {wj } . This 
grouping can be implemented with various different algorithms, e.g., averaging, 
self-organisation. Once we have constructed these sets we can use them to estimate 
a value of p(i) for data values recorded at t. 
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For data samples X(t), Y(t), Z(t) monitored at time t the vector u{t) = (X(t),Y(t) t Z(t)) 
is constructed. 

The distances D between u(J) and each element of {vi} and {wj} are calculated. 
This gives HJ distance calculations (as distance measure e.g., the Euclidean distance 
5 can be used). 

The smallest distance calculation result is noted. 

If this smallest distance was obtained when comparing u(f) with an element of {vi } 
then P(t) is assigned as being "high", if it was obtaining when comparing with an 
element of {wj} then P(t) is assigned as being "low". 

10 

It is also possible to use more than those two parameters described above. Accord- 
ingly, it is possible add measurement of another type cardiovascular activity in the 
first embodiment of the invention so that there are two signals and three parameter 
values for calculation of the probability index. In the second embodiment of the in- 

15 vention it is possible to derive a further parameter from the two signals so that there 
are three parameter values for calculation of the probability index. In the third em- 
bodiment of the invention it is possible to derive a further parameter from the two 
signals so that there are three parameter values for calculation of the probability in- 
dex, too. It shall be understood that it is also possible to perform further online 

20 measurements, whereupon three signals are acquired, and three or four or more pa- 
rameters are included in the calculation of the probability index. 

Example: 

In the following an example is described, which is based on an unpublished material 
25 concerning 23 patients, single-variable based estimators based on a PPG-signal, and 
a fEMG-signal for the level of analgesia could classify the patients 1 responses to 
skin incision with 65-74% correctness i.e. the physiological responses could dis- 
criminate the patients who moved or did not move as a response to skin incision. 
When the information in these variables is combined the correctness for detection of 
30 insufficient analgesia is improved to 83%. In the example of real-time data were 
recorded during abdominal hysterectomy surgery. After induction (Fig. 5D: Arrow 
A) with fentanyl 1 Jig/kg iv and propofol 1 mg/kg iv anesthesia was deepened with 
sevoflurane 8% in 100% oxygen via facial mask until endotracheal intubation (Fig. 
5D: Arrow B). Sevoflurane concentration was adjusted to equal 0.8 MAC (1.6% 
35 end-tidal). Surgery began 14 min after intubation with skin incision (Fig. 5D: Arrow 
C), after which fentanyl and propofol were administered. 



21 



For the end purpose at least two signals were continuously acquired from the pa- 
tients, in this example a pulseplethysmographic (PPG) signal and a facial electro- 
myographic (fEMG) signal. The instantaneous parameter values of R-to-R interval 
(= RRIpost and RRIpre), dicrotic notch height (= PtyNotchAmppost and PtyNotch- 
5 Amppre), fEMG power (= fEMGpost and fEMGpre) were continuously derived or 
calculated the heart rate and dicrotic notch being the used quantities. A nonlinear 
» combination of the above parameters according to the invention was calculated, as 
well as a linear combination of the same parameters as a comparison. The respec- 
tive instantaneous parameter values were measured real-time or online from each of 

10 the patients. The originally received raw signal, i.e. continuous curves of parameter 
values are not shown in the figures, but the curves normalized with a pre-incision 
reference level is visualized. Accordingly, Fig. 5 A discloses normalized R-to-R 
interval, Fig. 5B discloses normalized dicrotic notch height, Fig. 5C discloses 
normalized fEMG power, and Fig. 5D discloses index P for probability of patient. 

15 comfort obtained by mathematical nonlinear combination of the above mentioned 
normalized signals. In this mathematical function the ending -post means the latest 
parameter value and the ending -pre means the calculated historical data of the 
respective parameter values, which can be also called a reference value. As can be 
seen there used both subraction and division between the latest value and the 

20 reference value. The time axis is in seconds. 
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The actual non-linear equation for the mathematical index for probability was: 



P = 



-(-8,624^£^-2.672(/^o/cMm/>;^ 
+ g RRIpre fEMGpre 



where P is the probability of a movement as a response for noxius stimulus. This 
index is an indicator for potential nociception during anesthesia or sedation so that a 
high values of the index indicate high probability for nociception and low values in- 
dicate low probability of nociception. 

Signs for nociception and inadequate analgesia are clearly seen around intubation 
(Fig. 5D: arrow B) and incision (Fig. 5D: arrow C) in nonlinear combination 
parameter but not so clearly in individual variables. Note also recovery from 
anesthesia (Fig. 5D: arrow D) which is associated with an increase in the index. 



